Bioresource Technology 132 (2013) 423-426 



ELSEVIER 


Contents lists available at SciVerse ScienceDirect 


Bioresource Technology 


journal homepage: www.elsevier.com/locate/biortech 



Short Communication 

Sensitivity analysis of three-parallel-DAEM-reaction model 
for describing rice straw pyrolysis 

Junmeng Cai a,b *, Weixuan Wu a , Ronghou Liu a * 

a Biomass Energy Engineering Research Centre, School of Agriculture and Biology, Shanghai Jiao Tong University, 800 Dongchuan Road, Shanghai 200240, PR China 
b State Key Laboratory of Heavy Oil Processing, China University of Petroleum, 18 Fuxue Road, Beijing 102249, PR China 


HIGHLIGHTS 


► The three-parallel-DAEM-reaction model was successfully applied to describe rice straw pyrolysis. 
^ The kinetic parameters of the model were estimated by using the pattern search methods. 

^ The local sensitivity of the model to various parameters was analyzed in the first time. 

^ The parameters related to cellulose decomposition are more sensitive than other parameters. 
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The three-parallel-DAEM-reaction model was used to study the slow pyrolysis kinetics of rice straw 
based on thermogravimetric analysis (TGA) data. The kinetic parameters of the model were calculated 
using the pattern search method. A comparison between the predicted DTG data and experimental values 
showed good agreement. The influences of the kinetic parameters on the model for describing the 
experimental data of rice straw were analyzed by means of local parametric sensitivity analysis. The 
results indicated that the frequency factor and the mean value of the activation distribution for cellulose 
decomposition affect the model more strongly than other parameters, followed by the corresponding 
parameters for hemicellulose and lignin. The sensitivity of the model to the standard deviations of the 
activation energy distributions for all pseudocomponents is very slight. 

© 2012 Elsevier Ltd. All rights reserved. 


1. Introduction 

Biomass has great potential as a clean, renewable energy source 
that could dramatically improve our environment, economy and 
energy security. Therefore, the research on biomass energy conver¬ 
sion has attracted more and more attention in the world. Biomass 
can be converted to more valuable energy forms via a number of 
processes including thermal, biological, and mechanical processes 
(Bridgwater, 2012). Pyrolysis is a well-established thermal pro¬ 
cessing method for converting biomass into a liquid fuel that can 
substitute for fuel oil and make valuable chemicals (Bridgwater, 
2012 ). 

Biomass pyrolysis kinetics is required for the design and optimi¬ 
zation of pyrolytic reactors (White et al., 2011). Furthermore, bio¬ 
mass pyrolysis kinetics can produce submodels that can be coupled 
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with transport phenomena to describe the practical conversion 
process (Calonaci et al., 2010). 

Lignocellulosic biomass has a complex internal structure and 
contains three main components: hemicellulose, cellulose and 
lignin. Many models have been used to describe the pyrolysis 
kinetics of lignocellulosic materials (Bates and Ghoniem, 2012). 
The three-parallel-DAEM-reaction model is by far the most com¬ 
prehensive model for analyzing the pyrolysis kinetics of lignocellu¬ 
losic biomass (Varhegyi et al., 2011). The model assumes that 
lignocellulosic biomass contains three independently reacting 
pseudocomponents (hemicellulose, cellulose and lignin), and their 
pyrolysis kinetic behaviors can be described by the distributed 
activation energy model (DAEM). However, the sensitivity analysis 
of the model for describing biomass pyrolysis kinetics is still miss¬ 
ing. This is the scope of this work. 

2. Experimental section 

As a large agricultural country, China has huge potential of ligno¬ 
cellulosic biomass resource availability in the form of crop residues. 
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Rice is the second largest agricultural crop in China, and about 
200 million tons of rice straw are produced annually (Wang et al., 
2010). About 50% of rice straw can be used for bioenergy production 
according to the literature (Jiang et al., 2012). Therefore, rice straw, 
one representative type of lignocellulosic biomass, was used as the 
experimental material in this work. The rice straw samples were 
obtained from a farm located in Zibo, Shandong Province, China. 
The samples were used after screening to 100 mesh and dried for 
12 h at 378 K to remove moisture. Elemental composition and low¬ 
er heating value of rice straw are presented in Table 1. 

Slow pyrolysis of rice straw was performed using the thermo- 
gravimetric equipment Netzsch STA 449 C, which allows the mea¬ 
surement of mass changes and thermal effects. To achieve 
pyrolysis conditions, a nitrogen atmosphere was used. The flow 
rate of nitrogen was set to 25 mL min -1 for the purge gas and 
14mLmin -1 for the protective gas. Around 5 mg samples were 
spread on a Pt-Rh crucible. Thermogravimetric measurements 
were performed from room temperature to 900 K at the heating 
rate of 10 and 20 K min -1 . 



Fig. 1. Experimental and predicted DTG curves of rice straw pyrolysis at the heating 
rate of 10 K min -1 . 


3. Three-parallel-DAEM-reaction model applied to biomass 
pyrolysis kinetics 


In the three-parallel-DAEM-reaction model, it is assumed that 
lignocellulosic biomass contains three independently reacting 
pseudocomponents (hemicellulose, cellulose and lignin), and 
therefore, that the global thermal decomposition behavior reflects 
the individual behavior of these pseudocomponents, weighted by 
the composition. It is further assumed that the pyrolysis kinetics 
of each pseudocomponent can be described by the DAEM (Varhe- 
gyi et al., 2011). 

The equation of the model is given below: 
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In the above equation, a is the conversion degree, T is the tem¬ 
perature, p is the heating rate, E is the activation energy, R is the 
universal gas constant, k 0 is the frequency factor, E 0 is the mean va¬ 
lue of the activation energy distribution, o is the standard devia¬ 
tion of the activation energy distribution, c is the fraction of 
volatiles produced by a certain pseudocomponent, and the sub¬ 
scripts 1, 2 and 3 indicates the values related to hemicellulose, cel¬ 
lulose, and lignin, respectively. 

The unknown parameters of the model were determined by 
evaluating the experimental derivative thermogravimetric (DTG) 
data using least squares nonlinear methods. Therefore, the objec¬ 
tive function based on the DTG data was defined by the following 
equation: 
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where the subscript i refers to the data points used, nd is the num¬ 
ber of data points, (da/dT)exp, ij represents the experimental values 


Table 1 

Analytical characteristics of rice straw (dry ash free basis). 

C (% by wt) H (% by wt) N (% by wt) O a (% by wt) LHV (MJ kg -1 ) 

45.52 6.29 0.99 47.20 15.35 


at the heating rate j (j = 1 represents 10 K min -1 , and j = 2 repre¬ 
sents 20 K min -1 ), and (da/dT)cal, ij represents those calculated 
by Eq. (1) for a given set of parameters. 

There are many methods which have been used for solving the 
above optimization problem, such as the differential evolution 
algorithm (Santos et al., 2012), the direct search method (Gunes 
and Gunes, 2002), the simulated annealing method (Mani et al., 
2009), the flexible simplex method (Conesa et al., 1995), the Mul¬ 
tistart algorithm (Kirtania and Bhattacharya, 2012), and the pat¬ 
tern search method (Cai and Liu, 2008; Cai et al., 2011). The 
pattern search method is a derivative-free, direct search method 
and superior to other direct search methods in both robustness 
and number of function evaluations (Lewis et al., 2000). Therefore, 
the pattern search method was used to estimate the model param¬ 
eters. For this purpose, a program was developed in MATLAB for 
the implementation of the parameter estimation. 

The model parameters estimated from the above method were 
used to generate a DTG curve which was compared with the exper¬ 
imental one by the fit quality parameter. The fit between the 
experimental and estimated data can be found out by the fit qual¬ 
ity parameter stated below: 

FP(%) = 100 ( 3 ) 

E(d«/d T) mJ 
j=i 



Calculated by difference. 


Fig. 2. Experimental and predicted DTG curves of rice straw pyrolysis at the heating 
rate of 20 K min -1 . 
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Table 2 

Calculated kinetic parameters of rice straw pyrolysis, together with corresponding 
value of FP. 


Corn stover 

9 

Logtofcoj 

(fcoj/s- 1 ) 

Eq//kJ mol 1 

Gjl kj mol 1 

FP/% 

Primary hemicellulose 

0.13841 

13.134 

187.325 

4.142 

2.282 

Primary cellulose 

0.30193 

13.802 

204.393 

2.498 


Lignin 

0.22070 

15.984 

242.613 

38.782 



Table 3 

Sensitivity levels of parameters. 


Parameter 

Sensitivity level 

ko.i 

Medium 

£o,i 

Medium 

03 

Poor 

ko ,2 

High 

£o,2 

High 

02 

Poor 

ko,3 

Medium 

Eo,3 

Medium 

03 

Poor 


where (da/dT) mj is the maximum experimental DTG value at the 
heating rate j. Note that the lower the value of FP, the better is 
the fit (Trninic et al., 2012). 

Figs. 1 and 2 show the experimental DTG curves of rice straw 
pyrolysis at the heating rate of 10 and 20 K min -1 , respectively. 
From the two figures, it can be observed that the DTG curves con¬ 
sist of two visible peaks, which represents the decomposition of 
hemicellulose and cellulose, respectively (Varhegyi et al., 2011). 
Table 3 lists the calculated parameters of the three-parallel- 
DAEM-reaction model by means of the pattern search algorithm 
for rice straw pyrolysis. The predicted DTG curves along with the 
reaction rate curves of three pseudocomponents are also shown 
in the figures. The fit parameter, FP, is less than 2.3%, as shown 
in Table 2, indicating that the three-parallel-DAEM-reaction model 
fits the experimental data of rice straw pyrolysis very well. 

4. Sensitivity analysis of three-parallel-DAEM-reaction model 

Sensitivity analysis is the study of how the variation in the out¬ 
put of a model can be apportioned to different sources of variation 
in the model input and can provide a systematic procedure to esti¬ 
mate the accuracy and robustness of a model (Rao et al., 2009). 
Sensitivity analysis can be divided into two categories: global 
and local sensitivity analysis; the difference between them consists 





Fig. 3. Sensitivity of the three-parallel-DAEM-reactions model near the optimal kinetic parameters. 
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in the magnitude of the parameter changes to be studied (Sciaco- 
velli and Verda, 2012). A local analysis evaluates sensitivity relative 
to point estimates of parameter values while a global analysis ad¬ 
dresses sensitivity with regard to the entire parameter distribution 
(Confalonieri et al., 2010). In this work, we will focus on the local 
sensitivity analysis of the three-parallel-DAEM-reaction model 
for describing rice straw pyrolysis. 

The sensitivity analysis was obtained by varying the model 
parameter around the optimal value of the selected parameter 
and determining the corresponding change in the objective func¬ 
tion (2). To find the sensitivity of the model, the kinetic parameter 
is given ± 20% deviation from the optimal value. The objective 
function values were calculated and their dependencies on the ki¬ 
netic parameters are shown in Fig. 3. In these figures, the y-axis 
represents the relative objective function, calculated as the ratio 
of the objective function value at the deviated parameter to the 
objective function at the optimal parameter value. The x-axis rep¬ 
resents the ratio of the deviated parameter to the optimal param¬ 
eter value, which is called the relative parameter value. The 
relative objective function near one means that a parameter 
change has no significant effect on the dependent variable. On 
the other hand, as the relative objective function increases, the 
parameter changes in the dependent variable become more signif¬ 
icant (Santos et al., 2012). 

According to the results included in Fig. 3, the sensitivity levels 
of the kinetic parameters are depicted qualitatively in Table 3. In 
the table, “high”, “medium”, and “poor” represent the degree of 
sensitivity to the parameter, which is defined on the basis of the 
percentage variation of the objective function when a percentage 
change in a parameter is given around an optimal value. The fol¬ 
lowing inferences may be drawn from Fig. 3 and Table 3. 

(1) The sensitivity of the model to the standard deviations (cri, 
(72, cr 3 ) of the activation energy distributions for three 
pseudocomponents is slight. 

(2) The frequency factor (k 0>2 ) and the mean value (E 0i2 ) of the 
activation energy distribution for cellulose affect the model 
more strongly than the other kinetic parameters, followed 
by the frequency factors (k 0 ,i and k 0 , 3 ) and the mean values 
(£ 0> i and £ 0 , 3 ) of the activation energy distributions for hemi- 
cellulose and lignin. 

The sensitivity of the parameters (k 0<2 , £ 0 , 2 ) related to cellulose 
decomposition is very high, which is due to the fact that the 
decomposition range of cellulose is around the DTG peak of bio¬ 
mass pyrolysis. Some special consideration for the parameters 
k 0 ,2 and £ 0i 2 should be given in parameter estimation of the model. 

5. Conclusions 

(1) The pyrolysis kinetics of rice straw can be successfully 
described by the three-parallel-DAEM-reaction model. 

(2) The sensitivity of the model for describing rice straw pyroly¬ 
sis was investigated by local parametric sensitivity analysis. 
The results showed that the model is more strongly affected 
by the frequency factor and the mean value of the activation 
energy distribution for cellulose than other parameters, fol¬ 
lowed by the corresponding parameters for hemicellulose 


and lignin and the model has only a very slight sensitivity 
to the standard deviations of the activation energy distribu¬ 
tions for three pseudocomponents. 
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